Bose analogs of MIT bag model of hadrons in coherent precession 
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Recently it was demonstrated that magnon condensation in the trap exhibits the phenomenon of 
self-localization [J. When the number of magnons in the textural trap increases, they drastically 
modify the profile of the gap and highly increase its size. The trap gradually transforms from the 
initial harmonic one to the box with walls almost impenetrable for magnons. The resulting texture- 
free "cavity" filled by the magnon condensate wave function becomes the bosonic analog of the MIT 
bag, in which hadron is seen as a cavity surrounded by the QCD vacuum, in which the free quarks 
are confined in the ground or excited state. Here we consider the bosonic analog of the MIT bag 
with quarks on the ground and excited levels. 
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In Ref. [l| it was demonstrated that the Bose-Einstein 
condensation (BEC) of magnons in the flexible trap ex- 
hibits the phenomenon of self-localization with forma- 
tion of a box. This phenomenon is not unique in nature. 
Other examples of self-formation of a box-like trapping 
potential are the electron bubble in liquid helium and the 
MIT bag model of a hadron where the asymptoti- 
cally free quarks are confined within a cavity surrounded 
by the QCD vacuum. Here we discuss the analogy of 
magnon BEC with the MIT bag model. 

The MIT bag model has been used for construction of 
different hadrons, including mesons, baryons and even 
multiquark hadrons, such as tetraquarks Q and pen- 
taquarks Q. In the MIT bag model, free quarks are 
forced to move only inside a given spatial region, within 
which they occupy single-particle orbitals. MIT bag is 
described by the following energy whose minimization 
determines the equilibrium radius R of a given hadron: 



E(R)=J2 N *V m2 a c4 + 
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Here the first term is the kinetic energy of quarks with 
masses m a in the cavity of radius R, where parame- 
ters x a arc determined by the boundary conditions for 
fcrmions on the boundary of the bag. For the fermion in 
the ground state in the box, and in the ultra-relativistic 
limit of vanishing fcrmionic masses, m a — > 0, the param- 
eter x = 2.04. The second term is the potential energy, 
B is the so-called bag constant that reflects the bag pres- 
sure. At zero temperature, the bag constant B is the 
difference in the energy density between the false vac- 
uum inside the bag (the deconfinemcnt phase) and the 
true QCD vacuum outside (the confinement phase). In 
the non-relativistic limit, ignoring the term which does 
not depend on R, one gets 





FIG. 1: Original magneto-textural trap for magnons. The 
oscillator frequency lu~ of the well in the axial direction in 
(J5]) can be regulated by the field in the pinch coil. The fre- 
quency uj r in the radial direction can be adjusted by applying 
rotation, since the vortex-free superfluid flow or the array of 
rectilinear vortex lines created by rotation modifies the flare- 
out texture l(r). 



The same equation describes the electron bubble in su- 
perfluid 4 He, where m is the electron mass; x = n 2 
for the ground state level; the potential energy F(R) = 
(47r/3)i? 3 P + AitaR 2 , with P being the external pressure 
and a the surface tension. For extension to the multi- 
electron bubbles in superfluid 4 He see 

This consideration is applicable for magnon BEC. 
The original harmonic trap for magnons is schemati- 
cally shown in Fig. [T] jlj. The confinement potential 
f/|| (z) = r yH{z) in the axial direction is produced by local 
perturbation of magnetic field with a small pinch coil. In 
the radial direction the well U± (r) is formed by the cylin- 
drically symmetric flare-out texture of the orbital vector 
1 (direction of the orbital momentum of Cooper pairs). 
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It comes from the spin-orbit interaction energy 
U ± {r)\*\ 2 = ^L(l-l(r))\*\*, 

OWL 



(3) 



where ^ is the wave function of magnon condensate, and 
I = 1 • H = cos /3l describes the orientation of the orbital 
momentum 1 with respect to the direction H of magnetic 
held. On the side wall of the cylindrical container the 
orbital momentum 1 is normal to the wall, while in the 
center it is parallel to the axially oriented applied mag- 
netic field. This produces a minimum of the potential 
Ui_{r) on the cylinder axis. So the total confinement po- 
tential for the "f-field is 



U(r) = U {l {z) + U ± (r)=LJ L (z) + 
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P^l-lto). (4) 



Close to the axis the polar angle /3l of the 1- vector varies 
linearly with distance r from the axis. As a result, the 
potential U (r) reduces to that of a usual harmonic trap 
used for the confinement of dilute Bose gases Q : 



C/(r) = [/(0) + ^( Wz V +Wr V) 



(5) 



where we take into account that the axial trap is also 
close to harmonic. The harmonicity of the initial trap was 
checked by measurement of the spectrum of magnons - 
standing spin waves in the trap, which was found to have 
nearly equidistant levels. The equidistant levels of spin 
precession localised around the minimum of the external 
magnetic field have been derived in Ref. [|| using the full 
set of Lcggett equations for spin dynamics of 3 He-B. 

In the analogy with the MIT bag, the flexible texture of 
the orbital momentum 1 plays the role of either the pion 
field or the non-perturbativc gluonic field depending on 
the microscopic structure of the confinement phase. The 
trap is modified by pumped magnons due to spin-orbit 
interaction in © : the condensation of field ^ in the trap 
leads to the preferable orientation of 1 parallel to mag- 
netic field, I = 1, in the region of the trap. In other words 
the condensate repels the 1-field from the region, where 
magnons are localized. At large number Af of magnons 
in the trap the systems becomes similar to MIT bag with 
cavity free from the orbital field, which is occupied by 
magnons. So magnons, like quarks, dig a hole pushing 
the orbital field away due to the repulsive interaction, 
see Fig. [31 The main difference from the MIT bag model 
is that magnons are bosons and may macroscopically oc- 
cupy the same energy state in the trap, forming the Bosc- 
condensate, while in MIT bag the number of fermions on 
the same energy level is limited by the Pauli principle. 
The bosonic bag becomes equivalent to the fermionic bag 
in the limit of large number of quark flavors, when Af 3> 1 
quarks may occupy the same level. Without loosing gen- 
erality, we may consider the 2D approximation, i.e the 2D 
cylindrical trap. In the limit of large Af, the radius R of 



the cavity filled with Af magnons occupying the quantum 
state with radial number n r is determined by a balance 
of two terms in the total energy of the bag [l| : 



E(R,n r ) = Afe nr (R)+F(R) , e nr (R) 



h 2 X 2 



2m M R 2 



(6) 



The first term on the rhs is magnon zero-point energy 
in the cavity. It is the magnon number Af in the Bose 
condensate times the energy e„ r of a single magnon on the 
radial level n r in the cylindrical box with impenetrable 
walls. In NMR experiments only the energy levels with 
zero azimuthal quantum number are excited. They are 
measured as the shift of the frequency of the NMR peak, 
corresponding to excitation of a magnon, with respect to 
the Larmor frequency ujl'- 
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Here myi is as before the magnon mass, and the parame- 
ter x in (0) equals the n r +l-th root of the Bcsscl function, 
x = A nt +i, which corresponds to the proper boundary 
condition for the magnons populating the radial level n r 
in the impenetrable box. The potential energy F(R) in 
© corresponds to the pressure exerted to the bag by the 
field of the orbital texture, which is expelled from the 
bag. It is the difference in the energy of the orbital field 
texture with and without the cavity. 

Experimental results for the ground state magnon con- 
densate in [l| demonstrated that they can be reproduced 
by the phenomcnological equation ^ if one assumes that 
there is the scaling law F(R) oc R '". Minimization of the 
phenomcnological equation Eq. © with respect to R sug- 
gests that at large Af the radius of localization approaches 
the asymptote: 

R(Af) ~ a r (Af/Af c ) 1/{k+2} , AA»AA C , (8) 

where a r is the harmonic oscillator length in the orig- 
inal radial trap (at Af <C Af c ), Af c is the characteris- 
tic number at which the scaling starts, In experiments, 
the dependence of the transverse magnetization Ai± on 
the frequency shift Aui is measured. As distinct from 
the magnon number Af = JcPr^l 2 , where "J is the 
wave function of magnon condensate, the transverse mag- 
netization density represents the order parameter and 
is proportional to ty. The total magnetization is thus 
M± oc J d 2 r|*| oc Af 1 l 2 R. Since Aw oc 1/R 2 according 
to 0, one obtains A4j_ oc (Au)~ 1 ~ k / 4 . The measured 
transverse magnetization suggests that for large Af the 
scaling law is approached with k w 3. The magnetiza- 
tion estimated in numerical simulations also suggests that 
under experimental conditions the k = 3 scaling is the 
reasonable fit [l[ , which is just the scaling corresponding 
to the MIT bag model. 

In the above approach the energy consideration has 
been used, where the energy potential is obtained by av- 
eraging over fast precession. The numerical simulations 
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of the Q-ball, which used the full dynamical equations 
for spin and order parameter in one-dimensional textural 
trap, can be found in Q. 

Incidentally, for an atomic condensate in harmonic 
trap the radius R as a function of number of atoms 
at large Af also approaches the scaling in Eq. ([8]) with 
k = 3 [7J . This behavior results from the repulsive inter- 
particle interactions in the Thomas-Fermi limit. How- 
ever, this similarity in the scaling law k = 3 both with 
atomic condensate and with MIT bag model is accidental. 
Moreover, the k = 3 scaling is actually in disagreement 
with the typical bag models. If the main contribution 
to the pressure comes from the bulk, as it happens for 
the hadron model, then the energy F(R) should be pro- 
portional to the volume of the bag, and then for the two- 
dimensional radial trap one would expect F(R) = ttR 2 P 1 
i.e. k = 2. On the other hand, if the main contribution 
to the pressure comes from the surface tension, as it hap- 
pens for electron bubble in liquid helium at P = 0, then 
the energy F(R) should be proportional to the surface 
area, and for our 2D case this would give F(R) = 2it<tR, 
i.e. k = 1. The observed more soft behavior with approx- 
imate scaling law k w 3 in the 2D case reflects the flexi- 
bility of the orbital field, which is inhomogeneous outside 
the cavity. On general grounds, F(R) depends on sev- 
eral length scales: radius R of the bubble; radius R c of 
the cylindrical container, where the boundary conditions 
on the 1-texture are imposed; and the textural healing 
lengths: magnetic length (the thickness of the layer 
near the wall of the container in which the orientation of 
1 by magnetic field is restored) and the lengths related to 
the orientational effects of rotation and vortices on 1. 

We investigated numerically the scaling law emerging 
for the condensates formed on the ground and excited 
levels for textures in the rotating vessel of two different 
sizes. In the MIT bag model the shape of the bag is 
spherical if all quarks are in the ground state. When 
considering excited hadrons in which quarks occupy also 
the higher energy levels, the non-spherical shapes of the 
bag must also be considered. In our case this would cor- 
respond to the magnon bubble with broken axial sym- 
metry. However, so far in our NMR experiments we 
were able to excite and grow BEC condensates only with 
zero azimuthal quantum numbers, so that the shape of 
the bag with such condensates is automatically cylindri- 
cal. Figs. [5] and [3] demonstrate numerical simulations of 
the formation of the magnon condensates on correspond- 
ingly ground state level and on the excited radial level in 
the original harmonic trap formed by the vortex-free tex- 
ture in the experimental cell of radius R c = 3mm in the 
cryostat rotating as O = 0.8 rad/s. At large number of 
pumped magnons the original harmonic trap is modified: 
for all radial quantum numbers n r the trap gradually 
transforms to the cylindrical box free from the 1 texture. 

We checked how well the modified trap represents the 
cylindrical box with the wall impenetrable for magnons. 




radius, millimeters 



FIG. 2: Simulation of formation of the multi-magnon bubble 
with magnon condensate on the ground state level, n r — 0, 
in the cell of radius R c — 3 mm in the vessel rotating with 
angular velocity Q — 0.8 rad/s in the vortex-free state for dif- 
ferent number of pumped of magnons. (top): /3l is the polar 
angle of the orbital 1-vector, cos/3l(V) = l(r) ■ z = l(r). It 
plays the role of neutral field, which provides the trap po- 
tential for the charged field according to @. (bottom): 
Pm(t) is the tipping angle of precessing magnetization. It is 
related to the charged field - the wave function of magnon 
BEC - as ^ ~ sin(/?M/2) exp(iuit + ia), see Ref. Here uj 
and a are the frequency and the phase of coherent precession, 
which in the magnon BEC language are the chemical poten- 
tial and the phase of the Bose condensate correspondingly. 
The total transverse magnetization M± = J d V\H sin Pm (r) 
is connected with the number of magnons Af = J dV\^\ 2 oc 
f dV(l — cos/3m(^)). In figure, M denotes the normalized 
magnetization, M = V~ x f dV sin /3m (r). With increasing M 
and thus with increasing TV", the neutral field Pl(t) is more 
effectively repelled from the region of the magnon conden- 
sate due to repulsive interaction in ((3]). As a result, the har- 
monic trap gradually transforms to the box (top), while the 
wave function of the condensate gradually transforms to 
the Bessel function (bottom). 



For such a box filled with magnons on n-th level the ra- 
dius as a function of the frequency shift is given by (|T|). 
Figure |4] demonstrates the dependence of R on the fre- 
quency shift for different radial quantum numbers n r of 
the magnon condensate in the trap both from numerical 
simulations and from the phcnomenological equation ([7]) . 
One can see that with decreasing frequency the harmonic 
states of the condensate gradually approach the states in 
the box which obey Eq. ([7]). 

Both in experiments and in the numerical simulations 
in [f[, all the length scales were of the same order, and 
thus no really small parameter was available, which could 
justify the scaling law with k ~ 3. The true scaling 
behavior may only appear in some limit cases. One of 
the limit cases may take place in the large vessel rotat- 
ing with angular velocity f2 in a vortex-free state in the 
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FIG. 3: The same as in Fig. [2]but with the condensate formed 
at the first excited radial level in the trap, n r = 1. The fiat 
region with I = f develops at large transverse magnetization 
(top), indicating the formation of a box in which magnons 
are condensed on the first excited level (bottom). The wave 
function ^ of the condensate has a node. 
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FIG. 4: Log-Log plot of the radius of the magnon bubble 
as a function of the frequency shift for the magnon BEC in 
the ground state and on three excited levels. Solid lines - 
numerical simulations, dashed lines correspond to Eq. (J7J) 
describing the energy levels in a cylindrical box. 



regime, when R c 3> R S> ~> where £„ is the 
healing length related to counterflow |v s — v„| = Sir, 
which orients 1 in the plane. In the large vessel with 
counterflow the boundary of the vessel are not impor- 
tant, because the main orientational effect comes from 
the counterflow, so that within the cavity one has I = 1 
due to spin-orbit interaction of the 1-field with the 
condensate, while outside the cavity one has I = due to 
counterflow. As a result, the energy F(R) is simply the 
excess of the orientational energy of the counterflow in 
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FIG. 5: The same as in Fig. 0for the first excited state, but in 
the larger vessel (radius R c = 30 mm). On the plot only the 
central part of texture is shown: for r > 3 mm there is nothing 
interesting to be seen, since the texture practically reaches the 
constant value I = prescribed by the counterflow, while the 
wavefunctions are equal to zero. 

the cavity F(R) oc irR 2 (QR) 2 , i.e. the scaling law with 
k = 4 emerges in the limit of large vessel, or large coun- 
teflow. This gives for the radius of the magnon BEC bag 
at large Af the asymptote 

i?~A/" 1/6 , (9) 

We probed this regime in our simulations using the 
vortex free texture in the rotating vessel with large radius 
R c = 30 mm. Figure [5] demonstrates the texture and the 
condensate density in case of the magnon BEC on the 
first excited level. Qualitatively the behavior is similar 
to that in Fig. [31 obtained for the experimental cell with 
the radius R c = 3 mm. The power law in the dependence 
R(AT) of radius of the magnon BEC bag calculated for 
the ground state condensate, is 

R ~ TV" 0178 . (10) 

This power is more close to 1/6 w 0.167 in ([5]), which 
corresponds to k = 4, rather than to 0.2 for k = 3. 

The MIT bag model is used also for description of large 
(macroscopic) box with the dcconfincmcnt phase inside, 
which is filled with thermal quarks and gluons forming 
the quark-gluon plasma, and with the hadronic phase 
outside, which contains only massive particles, such as 
pions. Using macroscopic MIT box one can describe the 
confinement-deconfinement phase transition, see e.g.fioj]. 
In our case the analog of such transition occurs, when 
the magnon density reaches the value at which the ori- 
entational effect of the magnon BEC on the orbital field 
exceeds the competing orientational effect of the coun- 



5 



tcrflow in the whole vessel. This happens when 



sin/3 



M 



n 



L V ph 



(11) 



where v p b is the pair breaking velocity and £Il is the so- 
called Leggett frequency, which characterizes the spin- 
orbit interaction. When this threshold is reached, the 
texture (which is an analog of the pion field) is completely 
wiped out of the vessel, and the whole vessel becomes oc- 
cupied by the magnon BEC - the counterpart of plasma 
of free quarks and gluons. However, the coherent pre- 
cession in the whole vessel is known to experience Suhl 
instability - the catastrophic relaxation 11- 13j|. 

In conclusion, the extension of the MIT bag model 
to bosons allows us to describe the self-localization of 
magnon Bose-condensate. The magnon bag is the kind 
of the Q-ball [l4[ formed due to interaction between the 
magnon condensate described by the field \I/ with con- 
served charge Q = Af, and the orbital field y(r) = 1— l(r), 
which is the analogue of the neutral field [15J. In the 
process of self-localization the charged field modifies 
locally the neutral field x so that the potential well is 
formed, in which the charge M is condensed. Due to the 
repulsive interaction between these field in 3 He-B, in the 
limit of large Af the cavity is formed, which is void of 
neutral field and is filled with the charge field. As a re- 
sult one obtains the bosonic analog of the hadron in the 
MIT bag model: instead of the quarks in the ground state 
in a box potential, there are magnons which also fill the 
ground state. But since magnons are bosons, they form 
the Bose condensate in a box. We also applied the MIT 
bag model to the analogue of an excited hadron - the 
bag with magnons condensed in one of the excited levels 
in a box. 

This work is supported by the Academy of Finland and 
the EU - FP7 program (# 228464 Microkelvin) . 
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